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ABSTRACT 



We present VLBI and archival Karl G. Jansky Very Large Array (VLA) and 
Westerbork Synthesis Radio Telescope (WSRT) observations of the radio after- 
glow from the gamma-ray burst (GRB) of 2003 March 29 (GRB 030329) taken 
between 672 and 2032 days after the burst. The EVLA and WSRT data suggest 
a simple power law decay in the flux at 5 GHz, with no clear signature of any 
rebrightening from the counter jet. We report an unresolved source at day 2032 
of size 1.18 ± 0.13 mas, which we use in conjunction with the expansion rate of 
the burst to argue for the presence of a uniform, ISM-like circumburst medium. 

Subject headings: gamma rays: bursts 
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Introduction 



Gamma-ray bu rsts (GRBs) are sudden, short-lived (< 10 2 s) releases of energy of on 



the order 10 51 erg (IFrail et al. 



2001 



sudden, short 


-livec 


Bloom et all 


2003) 



20031 ) in a region of < 100 km. The 



central engine, thought to be a black hole with an accretion disk formed either by the 
collapse of a massive evolved star (collapsar) for long GRBs (> 2s) or by the merger of two 

re large, sudden 



1993|). 



stellar remnants (black holes or neutron stars) for short GRBs (< 2s). T 
energy release produces an e + e~,7 fireball in the form of collimated jets ( iPaczyriski 
Variability in the rate and velocity at which material is ejected from the central engine 
leads to internal shocks within the jet when fast-moving shells catch up t o slower ones. The 



resul ting 7-ray emission is responsible for the so-called prompt emission (IRees fc Meszaros 
I99J ) . 



Whereas prompt emission lasts for < a few minutes, GRB after glow emission ca n 



last far longer, up to years in t 



Berger et al 



2003 



re radio part of the spectru m, e.g. (IFrail et al. 



Taylor et al. 



2004 : 



2000; 



Pihlstrom et al 



20071 ). The long-lived afterglow 



emission is primarily synchrotron radiation is usually attributed to an in teraction with the 
shock ed external medium that exists behind the external forward shock |]Vleszaros fc Rees 
19971 ). More recently, perplexing features observed by Swift in the early X-ray afterglow 
have lent weight to a possible alternative scenario wherein a long-lived reverse shock 



decelerates slow ejecta at the back of the qrigina 



the shocked external medium (IGenet et al. 



2007 



outflow as it gradua 



Uhm fc Beloborodov 



ly ca tches up with 



2007V ). Whether the 



afterglow arises from the forward external shock or from a long lived reverse external shock, 
an external shock origin is strongly supported by measurements of the radio afterglow 
i mage size and its temporal evolution for the late time radio afterglow of GR B 030329 



rtTaylor et al. 



2004 



Oren et al. 



2004 ; 



Granot et al. 



2005 



Pihlstrom et al. 



20071). 



To explain the intensity and duration of observed GRB afterglow emission, the presence 
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of a circumburst medium must be taken into account. External shocks naturally arise 
when jet material interacts with the circum burst medium, convertin g kinetic energy of the 



fireball into particle energy and luminosity (IMeszaros fc Rees 



19931 ). These interactions 



are expected to be collisionless, mediated instead by the tangled and compressed magnetic 
fields that exist at the shock boundary. Electrons accelerating along magnetic field lines at 
the sho ck boundary produce the power la w emission spectrum characteristic of g amma-ray 



1998 



bursts (IMeszaros fc Rees 



Granot fc Sari 



1993 



Katz 



19941 ) . The spectral and temporal evolution fjSari et al. 



20021 ) of GRB afterglows, then, are governed by such factors as th e 



structure and dynamical evolution of the relativistic jet (for a review, see 



Granot 



20071) 



but also by the envi ronment (e.g., constant density versus a windlike density profile; 



Chevalier fc Li 



20001 ) . A detailed study of the size evolution of a GRB afterglow would 



make it possible to better constrain the density profile of the circumburst medium. 

Progenitors of long ga mma ray burst s are expected to be associated wit h the collapse of 



evolved stars (collapsars) ( jWoosley 



1993 



Paczynski 



1998 



Fryer et al 



19991 ). At least some 

GRBs should therefore exist in gas-rich environments with isotropic windlike p(r) oc r~ 2 
density profiles out to a few tenths of a par sec, where t here is a region of roughly uniform 



density corresponding to the stellar wind reverse shock (IRamirez-Ruiz et al. 



20051 ). Because 



the radio afterglow is detectable for months or years after the initial explosion, the jet would 
be expected to have sufficient time to traverse the stellar wind and enter the uniform shell 
of material that lies beyond. Several long GRBs, however, have shown light curves that are 
more characteristic of a jet propagating through a uniform p(r) = po medium for the entire 
observed duration of the radio afterglow. More direct and less model-dependent methods 
of determining the radial density profile could be used to better constrain the progenitor as 
well as to test the validity of the standard fireball model. 



In the simplest emission models, the flux density evolution of the GRB afterglow is 
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related to the density of the medium in which the GRB is p r esent 



afterg low from GRB 030329 was monitored by 



Taylor et al 



dWaxman 



fl2004h and 



1997). The 



Pihlstrom et al. 



( 120071 ) to obtain the first well-determined expansion rate for a GRB up to 247 and 806 days 
after the burst, respectively. The apparent diameter of the burst afterglow was shown to 
increase from 0.065 ± 0.022 mas to 0.172 ± 0.043 mas between day 25 and 83. Thereafter, 
the expansion rate decreased, with the burst reaching a size of 0.347 ± 0.090 mas at 
t = 806 days. The evo lution of the mean apparent expansion speed of the afterglow image 



( jPihlstrom et al 



20071 ) suggests a transition to non-relativistic expansion after about one 
year (see Fig. [3]). Models of the afterglow expansion predicted similar sizes at day 217 and 
day 806 in both the uniform density and the wind density profile case (see Fig. |4]), making 
it necessary to observe the afterglow at later times so that the preferred model for the GRB 
030329 circumburst environment could be unambiguously determined using this method. 
Here we present observations of the GRB 030329 afterglow taken up to 5.5 years after the 
initial burst in order to shed light on the validity of the different afterglow models. We 
also seek to place constraints on the density profile of the circumburst environment, the jet 
dynamics and structure, and the proper motion. 



GRB 030329 



The gamma-ray burst GRB 030329 was first detected by the Hi gh Energy Transie nt 



Exp 



orer 2 (HETE-2) satellite at 11:37 UTC on March 29th, 200 3 



20031 ). and was localized in the 



0.1685 ( jGreiner et al. 



optical bands by 



(TVanderspek et al. 



Peterson fc Pried ( 120031 ). Its redshift of 



20031 ). corresponding to an angular distance of 587 Mpc, makes 



it one of the nearest GRBs to EarthPJ The GRB 030329 radio afterglow is the brightest 



x The redshift of GRB 030329 is found to be z = 0.1685. Assuming a Acdm cosmology 
with H = 71 km s" 1 Mpc" 1 , tt M = 0.27, and fi A = 0.73, GRB 030329 is located at an 
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ever to have been observed, reaching a maximum flux density of 55 mJy one week after the 



burst at 43 GHz (IBerger et al 



2003ah . The relative proximity of GRB 030329 to the Earth 



and its high afterglow luminosity present an unprecedented opportunity to study the flux 



evolution of a GRB afterglow, as we 



resolved by the VLB A at 4.86 GHz flTavlor et al 



1 as the rate of expansion of the burst, which w as fully 



200^ 



2005 



P ihlstrom et al 



20071 ). 



3. Observations and Data Reduction 

3.1. Global VLBI 

On 21 October 2008, 2032 days after the burst, we used 9 antennas from the 
global VLBI Network (Effelsberg, Jodrell Bank Lowell, Medicina, Noto, Onsala, Torun, 
Westerbork, Green Bank, and Arecibo) to observe the 5 GHz continuum afterglow 
emission of GRB 030329. This array provided a baseline range between 266 and 6911 km, 
corresponding to fringe spacings between 1.8 and 47 mas. The total bandwidth observed 
was 128 MHz (8 IFs of 16 MHz each), dual polarization. The nearby source J1048+2115 
was used as the phase-reference calibrator with a 3 minute cycle (2 minutes on-source, 1 
minute on the calibrator). The source was tracked for 7.5 hours in total, which included 
3 hours with the Arecibo telescope. 

Data were recorded at 1 Gb/s rate, except for at the GBT telescope where the data 
recording rate was limited to 512 Mb/s. The data were correlated at the Joint Institute for 
VLBI in Europe (JIVE) with a 2 seconds averaging time and 16 channels, and the GBT 
data was correlated as 2-bit data with the magnitude bit set to high. This correlation 
scheme was successfully applied to our previous high sensitivity observations (Pihlstrom et 
al. 2007). 



angular distance of = 587 Mpc, with 1.0 mas corresponding to 2.85 pc. 
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Data reduction was done using the AIPS package. First the data were corrected for 
a slight position offset for the GBT antenna, and then for a position offset of the phase 
calibrator determined by the position of the ICRF source J1051+2119 observed every 30 



minutes. This ensures the position is measured re 



done in previous VLB1 observations of this source (jPihlstrom et al 



ative to the same po i nts as has be en 



2007 



Taylor et al. 



2005 



20041 ). A continuum map was produced using the inner 75% of channels in each IF, and 



resulted in a 6 //Jy beam 1 RMS noise level. 

A circular Gaussian was fitted to the measured visibilities and we derived an angular 
diameter of 1.18 ± 0.13 mas. The error in this value comes from uncertainties in the size of 
the gaussian fit, and arises due to the RMS noise contained in the image. It should not be 
confused with the size of the source, which is not resolved. Rather, it is a measurement of 
the angular resolution of the image, and can only serve as an upper limit to the source size. 
The position obtained for GRB030329 at day 2032 is R.A. = 10 h 44 m 49 s .95944 and Dec = 
21°31'17.4372". 



3.2. VLA and WSRT Archival Observations 



Here we present late-period 5 GHz observations, taken from the NRAO and WSRT 
data archives, that were made between 672 and 1773 days after the burst (see Tabled]). 
The VLA data were reduced in the AIPS software package. The bright calibrator 3C286 
was used for absolute flux calibration, and J1051+213, which was separated from the 
source by 1.6°, was us ed to determine the c omple x an tenna gains. The W S RT d ata were 



originally published in 



van der Horst et al. 



(120051) and 



van der Horst et al 



(|2008|). WSRT 



data appearing in those two works covering the same time span as the VLA data that we 
have taken from the NRAO archive were re-analysed to ensure uniformity of analysis. The 
re-analysed data were reduced in AIPS, with 3C286 serving as the primary calibrator and 
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Date 


At 


Frequency 


Flux Density 


Instrument 




(days) 


(GHz) 


OJy) 




2005 Jan 29 


671.90 


4 A 


% 


324 ± 28 


VLA-B 1 


2005 Mar 31 


732.74 


U 


S6 


196 ±31 


VLA-A 2 


2005 Apr 07 


739.73 


U 


S6 


270 ± 29 


VLA-C 1 


2005 May 14 


777.06 


4 A 


55 


276 ± 38 


WSRT 5 


2005 Jun 08 


801.53 


U 


S6 


279 ± 27 


VLA-A 2 


2005 Nov 27 


974.51 


U 


SO 


197 ±40 


WSRT 6 


2006 Mar 22 


1088.7 


u 


S6 


169 ±31 


VLA-A 3 


2006 Apr 30 


1128.0 


u 


SO 


149 ± 34 


WSRT 7 


2008 Feb 04 


1773.3 


i.l 


S6 


70 ± 17 


VLA-A 4 



Table 1: VLA and WSRT observations of the GRB 030329 radio afterglow at 5 GHz. The 
data were taken under VLA project codes AF414 1 , AK583 2 , AS864 3 , and TYP100 4 , and 
WSRT sequence numbers 10502389 5 , 10506025 6 , and 106 2115 7 . W SRT data appearing 



i n thi s table correspond to data from 



van der Horst et al. 



(12005h and 



van der Horst et al. 



( 120081 ) that were re-analysed to ensure uniformity of analysis across all data points. 



either 3C48 or 3C147 serving as the secondary calibrator. Flux densities derived from the 
WSRT data that were re-analysed in this work appear in Table [H 



Using the VLA data in c 



VLA data ( 



Frail et a 



( Ivan der Horst et al 



2005 



2005 



uded in this work i n conjunction with pr eviously-published 



Berger et al 



2003a 



van der Horst et al 



Pihlstrom et al 



2008 



20071 ) and WSRT data 



we find that the decrease in flux 



density from day 59 to day 1773 (Fig. [T]) is well described by a power law F v oc t a with 



the tempora 
obtained by 



index a = 1.27 ± 0.03. This is consistent with the value of a = 1.23 ± 0.03 



Pihlstrom et al. 



(120071 ) for the time range of 59 to 806 days after the burst. 
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Fig. 1. — The GRB 030329 5GHz light curve can be represented as a decaying power law with 



a = 1.27 ± 0.03 after day 59. Red circ 
square s are data taken from this work, 



es are data taken fr o m the VL A archive, while blue 



van der Horst et al 



d2005|), and 



van der Horst et al. 



(120081 ). The vertical lines at the bottom of the figure denote dates at which the afterglow 
size was measured with VLBI. The green cross represents the VLBI data point presented in 
this work. 

4. Results 



4.1. Flux Density Evolution 

Our archival VLA and WSRT observations show that the evolution of the afterglow 
flux density follows a single power law drop-off of the form F v oc t~ a with a temporal 



index a = 1.27 ± 0.03 (see Fig.g ) 



previous work by 



Pihlstrom et al. 



rom day 59 through day 1773. This is consistent with 
(120071 ). who found that the drop-off was F v oc r L23±0 - 03 



between days 59 and 621. Our VLBI measurement on day 2032 is also consistent with 
the overall power law, with an observed flux density of 76 ± 10 /iJy closely matching the 
predicted flux of 73 /iJy. VLBI instruments, lacking short baselines, are much less sensitive 
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to extended emission than smaller interferometers such as the VLA and the WSRT. In 
order for measurements of the flux density made with the VLA and the WSRT to agree 
with measurements taken from VLBI observations, there cannot be significant extended 
emission. 

The electrons that are injected into a gamma-ray burst shock are generally assumed to 



follow a power law distribution in energy iV(7 e ) oc 7 e p . The temporal index can 



De used to 



20081 ). In the 



determine the value of the electron power law index p (Ivan der Horst et al. 
non-relativistic regime, the temporal index and the electron spectral index, p, are related 
by a = 3(5p — 7)/10 in the homogenous case and by a = (7p — 5)/6 in the wind case. Our 
value of a = 1.27 ± 0.03 yields p = 2.25 ± 0.02 for a uniform medium and p = 1.80 ± 0.03 for 
a wind. If the circumburst density is assumed to follow a power law in radius p cxt oc r~ k , 
then the slope power law k can be calculated directly using 



5a -15/3 + 3 
a -4/3 + 2 ' 



wher e is the spectral in dex 



from 



van der Horst et al 



van der Horst et al 



( 12008 



2008 1. Using the value of = 0.54 ± 0.02 



one finds k — 1.1 ± 0.2. It is important to note, however, 
that the value that is obtained for p (and, therefore, k) using this method is extremely 
sensitive to the time range that is used in the calculation. The transition of the jet to 
non-relativistic expansion leads to a steeper decay phase than in the Newtonia n regime, and 



will lead to a higher value of p. It is contended by 



van der Horst et al. 



(120081 ) that the jet 



becomes non-relativistic at t ~ 100 days, and so any data before this time is excluded. With 
this choice of cutoff time, they obtain a value of p — 2.12 ± 0.02 and k = 0.331q4i. They 
point to the fact that they obtain larger values of p if they include earlier data to support 
their choice of the cutoff time. Setting the cutoff to t = 312 days, however, yields an even 
larger value of the electron power law index of p = 2.43 ± 0.07 (k = 1.9 ± 0.19), even though 
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the burst is almost certainly entirely non-relativistic by this time. Setting the cutoff time 
to t = 426days, on the o ther h and, yields p = 2.14 ± 0.07 (k = 0.46 ± 0.44. It is also noted 



by 



van der Horst et al. 



(j2008l ) that the value of x 2 decreases significantly if data taken 
prior to t = 100 days are excluded, which they use as additional support of their choice 
of cutoff time. Choosing a cutoff time of t = 312 days, however, will yield an even lower 
value of x 2 tha n a cutoff time of t = 100 d ays. We find that the value of k = 1.27 ± 0.03 is 



consistent with 



van der Horst et al 



(120081 ) given the large uncertainty in k that is obtained 



by calculatin g it directly from the 



however that 



van der Horst et al 



temp oral and spectral slopes. It should also be noted 



(120081 ) used multi-frequency data to obtain their value of 



p, meaning that they had more independent data points with which to calculate it. 



4.2. Size and Expansion Rate 



The entire history of expansion for GRB 030329 is shown in Fig. [2J The first 
measurement at 15 days comes 



scintillation (IBerger et al. 



*rom a model-dependent estimate of the quenching of the 
2003a ). The uncertainties on this size estimate are large due to 
the dependence of the measurement on estimated properties of the interstellar medium 
along our line-of-sight. A 2a upper limit on the afterglow image size of 1.31 mas, or 
1.2 x 10 19 cm, was obtained for day 2032. 

The apparent expansion speed is defined as /3 app c The decline in (/3 app ) with time is 
shown in Fig. [3l where (/3 app ) is found using 



</3app) 



(l + z)R ± 9 R d M 



ct 



ct 



(2) 



where 6r = R±/d,A = (1 + z)R±/cIm and R± are the angular and physical radius of the 
image, respectively, du is the proper distance to the source, z is the source's cosmological 
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redshift, t is the time of o bservation, and c is the s peed of light. Previous work by 



Pihlstrom et al 



f l2007h and 



Taylor et al 



( 12004 



20051 ) shows a continuous decrease in the 



apparent expansion rate of the image up to day 806, likely due to a deceleration of the 
afterglow shock due to interaction with the circumburst medium. The mean afterglow 
apparent expansion rate up to day 806 was (/? a pp) — 0.9 ± 0.2. Our upper limit of 1.31 
mas at day 2032 corresponds to a mean apparent expansion rate of ((3 app ) < 1.3. It should 
be noted that the intrinsic image surface brightness profile of the afterglow is expected 



to have a small effect on the va 



image may be limb-brightened (IGranot et al 



ue obtained 



or the diameter of the after glow 2R±. The 



1999 



Granot fc Loeb 



20011), which could 



suggest a possible correction factor to the values used in this paper of ~ 1.4. For an 
in- depth discussion of the effects of surface brightness profiles on the afterglow diameter, 



sec 



Pihlstrom et al. 



(I2007f ). 



4.3. Proper Motion 

Solving for proper motion using all the VLBI observations to date, we derive 
/^r. a. = 0.0044 ± 0.054 mas yr" 1 and //dec. = —0.0031 ± 0.033 mas yr _1 , or an angular 
displacement over 2032 days of 0.030 ± 0.35 mas (Fig. [5D. These observations are consistent 



with those reported by 



Taylor et al 



( 12004 



20051 ) and 



Pihlstrom et al 



(120071 ). and impose 



a stronger limit on the proper motion. The implied 2a limit on the proper motion in the 
plane of the sky is 0.067 mas yr" 1 (corresponding to 1.1 pc) in 2032(1 + z)^ 1 rs 1739 days, 
or 0.73c (la). 



In the relativistic fireb all model, a shift in the flux centroid is expected owing to the 



spreading of the jet eiecta (ISari 



substantial ( IGranot fc Loeb 



19991 ). For a jet viewed off the main axis, the shift can be 



20031 ). However, since gamma rays were detected from GRB 



030329 it is likely that we are viewing the jet largely on-axis. 
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Fig. 2. — The apparent expansion of GRB 030329 derived from measurements and limits on 



the an gular size as a function of time. The la upper limit at day 51 is from 



Taylor et al. 



(12004). as are the measurements on days 25 and 83. The mea surement on day 217 i s from 



Taylor et al 



( 120051 ). The measurement on day 806 comes from 



Pihlstrom et al 



(120071 ). The 



measurement on day 15 comes from a mo del- dependent estimate based on the quenching of 



the scintillation (IBerger et al. 
this work. 



2003aJ). Finally, the la upper limit on day 2032 comes from 



The angular shift of the flux centroid is approximately q8 R/D ( jSari 



19991 ). where q is 



a measurement of the angular separation between the observer and the center of the jet in 
units of the half-opening angle of the jet, 6q is the initial half-opening angle of the jet, R 
is the radius of the emitting region, and D is the angular distance to the observer. The 
maximum angular shift is found by setting q = 1, in which case the observer is at the edge 
of the jet. For a typical value of the initial jet half-opening angle of 6q = 0.2 mas, our 
upper limit on the emitting region size of R < 1.2 x 10 19 cm, and a luminosity distance of 
D = 587 Mpc, the maximum angular shift in the flux centroid is 0.16 mas. Although this 
estimate is larger than our upper limit of 0.067 mas, it can easily be made consistent by 
either decreasing the value of 0q or of q. Indeed, the detection of a gamma-ray afterglow 
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makes it likely that we are positioned near the center of the jet and q <^ 1. 



Proper motion in the cannonball model originates from the su per luminal motion 



of plasmoids ejec t ed du ring a supernova explosion with r ~ 1000 (IDado et al. 



2003|). 



Par &: de Rujulal ( 120031 1 predicted a displacement of 2 mas by day 80 assuming plasmoids 
propagating in a constant density medium. This estimate was revised downward to 0.55 
mas by incorporating plasmoid interacti ons with density inhomogeneities at a distance of 



~ 100 pc within a wind-blown medium (IDado et al 
consistent with our proper-motion limits. 



20041 ) . Neither variant of this model is 



5. Discussion 
5.1. Density Profiles 

Fig. H] demonstrates how the measured evolution of the image size for GRB 030329 



Granot et al. 



compa res with the predictions of a suite of theoretical models developed by 
( 120051 ). The data constrain the external density profile at radii R > 10 18 cm. Our 



observation at day 2032 suffered from a very low signal to noise, resulting in a large 
uncertainty. Additionally, the source is unresolved, meaning that we can only report an 
upper limit of 1.31 mas (2a). 

Looking at Fig. HI it can be seen that the upper limit that we obtain for the image 
size on day 2032 cannot be used by itself to distinguish between a uniform and a windlike 
medium. Another, tighter constraint on the angular size can be inferred from the apparent 
afterglow expansion rate (/3 app ) (Fig. [3]) on physical grounds. The afterglow expansion rate 
at day 806 is (/3 apP ) = 0.9 ± 0.2. Although it is plausible that (/3 apP ) might have remained 
constant between days 806 and 2032, there is certainly no plausible physical mechanism 
that would cause it to increase. The maximum afterglow size at day 2032 can therefore 
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Fig. 3. — Evolution of the average apparent expansion velocity derived from direct size 
measurements, and ass uming a Gaussian intrinsic surface brightness profile. The la upper 



limit at day 51 is from 



Taylor et al 



meas urement on day 217 i s from 



2001 



, as are the measurements on days 25 and 83. The 



from 



Pihlstrom et al. 



Taylor et al. 



2005 



The measurement on day 806 comes 



20071 . The measurement on da y 15 comes 



rom a model-dependent 



2003al ). Finally, the la 



estimate based on the quenching of the scintillation (IBerger et al. 
upper limit on day 2032 comes from this work. 

be calculated by assuming a constant expansion rate of (Papp) = 0.9 over this period. The 
result is a maximum afterglow size of 0.90 mas (8.0 x 10 18 cm), significantly lower than the 
measured upper limit of 1.31 mas (1.2 x 10 19 cm) and at the low end of the range of values 
that are compatible with the wind models. 

Because we only have an upper limit on (/3 a pp) a t day 2032, we cannot immediately 
rule out a large deceleration due to the transition of the jet into a more dense environment 
such as the ISM that lies beyond a star's wind termination shock. A non-relativistic jet 



the light curve (IMesler et al 



2012 



Ramirez- Ruiz et al. 


2001. 


2005; 


Dai & Lu 


2002) 
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jet transitions to the non-relativistic regime when (/? a pp) — 2 (jPihlstrom et al.l 120071 ). or 
after approximately 1 year. The lack of any rebrightening in the 5 GHz light curve argues 
against an increasing density. The observations of the afterglow size at days 217 and 806 
are consistent with a uniform medium. With limits on the possible afterglow expansion 
rate making a windlike medium unlikely and a lack of any rebrightening in the light curve 
making a transition to either a denser uniform medium or a k > density profile unlikely as 
well, it is most probable that the jet is propagating through a continuous, uniform medium. 
n a follow-up paper we pla n to use state-of-the-art relativistic hydrodynamic simulations 



De Colle et al 



2012 



20111 ). which would provide a significantly more realistic modeling of 



the GRB jet dynamics. 



5.2. Rebrightening 



Radio rebrightening of GRB 030329 has been predicted by 



Granot fc Loebl (J2003J) as 



the counter jet becom es non-rela t ivistic and its emission is no longer strongly beamed away 



from us. At 15 GHz, iLi fc Song! (|2004| ) predicted that the flux density of the afterglow 
would be 0.6 mJy 1.7 yr after the burst. The actual 4.86 GHz flux density was ~ 0.33 mJy 
after 1.7 yr (see Fig. [T]). Using the relationship F v oc i/~ 0,6 , which is valid between these 
two frequencies, one would expect to find a corresponding flux density at 15 GHz of ~ 0.17 
mJy, a factor of ~ 3.3 lower than the prediction. A more recent model by 



van Eerten et al 



(120101 ) predicts that the contribution to the flux by the counter-jet should be significantly 
smaller and should occur much later at ~ 600 days. We find that the late time 4.9 GHz 
light curve, from 59 days to 2032 days, is consistent with a single power law decay of 
F v oc £- L27±0 - 03 (see Fig. [I]), which implies that, up to ~ 2032 days after the burst, there is 
no significant contribution to the observed flux density from the counter jet. Consequently, 
we find no clear evidence for a rebrightening up to 5.5 years after the burst. We note that 
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detailed numerical simulations show that the rebrightening due to the counter-jet beco mes 



De Colle et al. 



20121 ). In 



much less prominent for a wind-like density profile (see Fig. 13 of 
such a case it might still be consistent with the observed single power-law flux decay in the 
late time radio afterglow light curve of GRB 030329. 



6. Conclusions 

We present measurements of the 5 GHz flux density and image size of the GRB 
030329 radio afterglow taken over a period of 5.5 years. These observations clearly 
demonstrate that the expansion rate has decreased over time, with a transition to the 
non-relativistic regime at ~ 1 yr. After approximately day 59, the afterglow flux density 
follows a power law of F v oc t~ a , with tem poral index a = 1.27 ± 0.03, which agrees with 



the value of a = 1.23 ± 



van der Horst et al. 



03 obtained by (IPihlstrom et al 



20071 ). Using the method of 



(120081 ). the electron power law index p is found to be p = 2.24 ± 0.02 
for a uniform medium, which does not agree with their value ofp = 2.12±0.02. The value 
determined for the electron power law index is found to be highly sensitive to the time 
range used to calculate it, making the temporal slope only capable of providing a quick 
estimate of the electron power law index. 

A rebrightening of the source was expected as the counter jet became non-relativistic, 
however, no rebrightening was detected up to 5.5 years after the burst. Numerical 
simulations suggest that any rebrightening in a windlike medium would be more difficult to 
detect as compared to a uniform medium. However, a windlike medium is not favored by 
the measured evolution of the afterglow image size. A possibility also exists for asymmetry 
between the jet and the counter jet. If the counter jet had a lower initial energy or 
encountered a less dense external medium, then it would produce a smaller rebrightening 
than expected. 



- 18 - 



An upper limit of 0.067 mas yr _1 is found for the proper motion. Consequently, the 
proper motion of the flux centroid is constrained to be smaller than the diameter of the 
image, which is consistent with the fireball model but not with the cannonball model. 

The upper limit that we obtained for the source size at day 2032 (Fig. H]) is, by itself, 
insufficient to conclusively determine the nature of the circumburst medium through which 
the jet is propagating. We argue that, on physical grounds, the mean expansion rate (/3 app ) 
cannot increase between two epochs. The value of (Aipp) at day 806 was 0.9 ± 0.2, which 
corresponds to a maximum angular size at day 2032 of 0.90 mas, which is not consistent 
with a windlike medium. We therefore argue that the expansion rate of the burst favors a 
uniform medium 
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Fig. 4. — Tentativ e fits of theoretical models for the evolution of the source size (from 



Granot et al. 



20051 ) to the observed image size (of diameter 2R \) of the radio afterg low of 



GRB 030329 up to 83 days. The la upper limit at day 51 is from 



Taylor et al. 



the measurements on days 25 and 83. The measurement on day 2 17 is from 



2004 as are 



and the measurement at 806 days is from 



Pihlstrom et al 



Taylor et al. 



2005 



20071 ). The 2a upper limit at day 



2032 comes from this work. All direct measurements are plotted with la error bars. In model 

1 there is relativistic lateral spreading of the GRB jet in its local rest frame, while in model 

2 there is no significant lateral expansion until the jet becomes non-relativistic. The external 
density is taken to be a power law with the distance r from the source, p ext = Ar~ k , where k 
= for a uniform external density while k = 2 is expected for a stellar wind environment. 



-24- 



CD 
CO 



03 



1.0 



0.5 



0.0 



o 

CD -0.5 



• i -r 



80ct21 ' 




03Apr6 



'03Apr1 



0.5 0.0 -0.5 

R.A. Offset (mas) 



-1.0 



Fig. 5. — The positions derived from the observations in eight epochs relative to the first 
determination on April 1st at 8.4 GHz. Observations at multiple frequencies at a given 
epoch have been plotted separately since they are independent measurements. A circle with 
a radius of 0.29 mas (2a) is shown to encompass all measurements except for the one taken 
on 08 Oct. 21 (day 2032) which lies just outside the circle's edge, and those taken w ithin the 



2004h . Taken 



first eight days at 5 GHz, which suffer from systematic errors (ITaylor et al. 
together these observations provide a constraint on the proper motion of < 0.067 mas yr _1 
over 2032 days. 



